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EXECUTIVE  SUMMARY 


Detailed  microstructural  alterations  (damage)  in  Indiana  limestone  and  one  type  of 
concrete  under  different  loading  conditions  have  been  observed  using  liquid  metal  alloy 
preservation  techniques.  The  technique  of  using  low  melting  point  alloy  is  proved  to  be 
effective,  convenient  and  reliable.  The  basic  observations  show  the  following: 

The  basic  mechanism  for  initiation  of  extensile  microcracks  in  a  material  under  compres¬ 
sion  are  the  localized  tensile  stresses  induced  by  heterogeneities.  These  localized  tensile 
stresses  include  bending  of  long  grains,  Herzian  contact  tension  between  adjacent  grains, 
mismatching  elastic  properties,  boundary  discontinuities  and  pore  space  and  soft  inclusion 
induced  tension.  Once  the  microcracks  are  generated,  regardless  of  the  initial  orientation, 
they  generally  propagate  in  the  direction  of  the  maximum  compression.  The  propagation 
of  microcracks  is  controlled  by  the  combination  of  local  tensile  stress  and  the  overall 
confining  stress,  before  any  interaction  takes  place.  It  is  hypothesized  that  the  generation 
and  propagation  of  microcracks  also  change  the  local  stress  distribution,  which  in  turn 
affects  the  generation  of  other  microcracks.  Preliminary  statistical  results  for  both  the 
concrete  and  limestone  show  that  the  average  length  of  the  microcracks  decreases  with  an 
increase  in  confining  stress.  Average  aperture  of  the  microcracks  also  decreases  with  an 
increase  in  confining  stress.  The  aspect  ratio  (the  ratio  between  the  aperture  of  a  crack 
and  its  length)  remains  relatively  constant  for  concrete  and  fluctuates  for  Indiana  limestone. 
Most  of  the  microcracks  are  parallel  or  subparallel  to  the  direction  of  the  maximum 
compressive  stress.  The  average  orientations  of  microcracks  in  all  of  the  samples  tested  are 
within  15  degrees  of  the  maximum  compressive  direction.  The  density  of  microcracks  (the 
number  of  microcracks  per  unit  observed  area)  increases  with  an  increase  in  confining 
stress.  The  relationship  between  microcrack  distribution  and  applied  stresses  could  not  be 
quantified  in  this  phase  of  the  research.  The  macroscopic  failure  of  Indiana  limestone 
samples  under  unconfined  conditions  is  usually  in  axial  splitting,  and  under  high  confining 
pressure  in  the  form  of  shear  bands.  The  axial  splitting  of  specimens  is  usually  the  result 
of  one  or  several  extensile  cracks  propagating  all  the  way  to  the  ends  of  the  samples.  The 
shear  bands  are  observed  to  be  formed  from  many  short  extensile  microcracks  which  are 
subparallel  to  the  direction  of  the  maximum  compression  and  with  their  centers  slightly 
offset  from  each  other. 

Ultrasonic  wave  velocities  were  measured  concurrently  during  static  compression  tests. 
Artificial  “dynamic”  stress-strain  curves  were  reconstructed  from  the  dynamic  modulus 
calculated  from  the  velocities.  “Dynamic”  stress-strain  curves  reconstructed  from  ultrasonic 
velocities  showed  similar  behavior  to  the  concurrently  generated  static  stress-strain  curves. 
Both  static  and  dynamic  moduli  decreased  near  peak  stress  (strength)  level,  indicating  the 
generation  and  interaction  of  microcracks  (damage).  This  suggests  that  if  a  relationship 
between  the  “dynamic”  and  static  stress-strain  curves  can  be  established,  a  non-destructive 
method  such  as  ultrasonic  detection  can  be  used  to  quantitatively  determine  and  detect  the 
static  behavior  of  a  geomaterial  and  the  damage  level. 


Using  energy  conservation  principles,  theoretical  derivations  of  effective  elastic  properties 
of  an  isotropic  material  containing  axisymmetrically  distributed  microcracks,  showed  that 
moduli  of  a  material  decrease  with  an  increase  in  damage  density.  This  is  verified  by 
independent  concurrent  triaxial  compression  tests,  ultrasonic  velocity  measurements  and 
microstructural  observations.  These  results  show  that  if  a  one-to-one  relationship  can  be 
established  between  the  theoretically  derived  formulations  and  microcrack  distribution,  it 
could  be  possible  to  extrapolate  microcrack  density  at  any  stage  of  loading.  This  requires 
a  calibration  of  the  theoretical  formulations  given  the  stress  conditions  and  quantification 
of  microcrack  distribution  under  the  same  conditions.  In  a  Phase  II  proposal  for  this 
research,  submitted  under  separate  cover,  this  calibration  and  microcrack  density 
extrapolation  is  proposed. 
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1  Administrative  Background: 

The  following  reports  on  a  Phase  I  SBIR,  awarded  to  Terra  Tek,  Inc.  The  specifics 

are: 

Topic  Number  AF  90-186 

Title:  Statistical  Study  of  Stress-Induced  Microstructures  in 

Selected  Rock  and  Concrete 

Submitted  by:  Terra  Tek,  Inc. 

420  Wakara  Way 

Salt  Lake  City,  Utah  84108 

Submitted  to:  AFOSR/XOT 

Building  410,  Rm  A-119 
Bolling  AFB,  DC  20332-6448 

2  Technical  Background: 

An  abbreviated  summary  of  the  work  performed  is  as  follows; 

2.1  Issues  to  be  Addressed 

In  this  proposed  study,  the  following  questions  will  be  addressed: 

•  What  are  the  mechanisms  for  the  generation  of  stress  induced  microcracks  in 
granular  and  porous  brittle  material  such  as  limestone,  sandstone  and  concrete,  when  these 
materials  are  subjected  to  overall  compression?  What  are  the  effects  of  inher¬ 
ent/pre-existing  microstructures  in  these  porous  material  on  the  generation  of  microcracks 
under  loading? 

•  What  are  the  major  factors  controlling  the  propagation  of  fractures  once  they  are 
initiated?  Are  they  related  to  the  microstructures  or  stress  conditions  or  both?  Can  these 
effects  be  quantified? 

•  What  are  the  interactions  between  microcracks?  How  do  local  stress  concentrations 
around  the  cracks  affect  the  generation  and  propagation  of  other  cracks? 

What  is  the  relationship  between  microcrack  density,  length  and  aperture  and  stress 
conditions? 
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•  What  are  the  mechanisms  of  macroscopic  failure  (axial  splitting,  shear  faulting  and 
drumming  out)  in  test  specimens?  How  do  microcracks  interact  to  cause  these  macroscopic 
failures? 

•  How  is  the  dilation  (expansion  in  the  minimum  stress  directions)  of  specimens 
related  to  the  generation  and  propagation  of  microcracks?  What  role  does  pore  collapse 
play? 

•  What  are  the  effective  elastic  moduli  when  microcracks  are  generated,  especially  the 
anisotropic  elastic  moduli  in  relation  to  the  microcrack  density,  length  and  aperture?  What 
are  the  anisotropic  Poisson’s  ratios? 


2.2  TASK  1:  Experimentation,  Microstructural  Observation  and 

Statistical  Analysis 

Objective:  To  statistically  study  and  quantify  the  generation,  propagation  and  interaction 
between  microstructures  generated  under  different  stress  conditions  and  the  mechanisms 
of  generation  of  these  stress-induced  microstructures. 

Technical  Approach:  Experiments  were  performed  on  one  type  of  rock  and  one  type  of 
concrete.  Specimens  were  subjected  to  uniaxial  and  triaxial  compression,  and  tension.  The 
microstructures  before  and  after  the  tests  were  examined  using  electron  microscopy. 

•  Rock  and  concrete  specimens,  one  inch  in  diameter  and  two  inches  long,  were 
endground  flat  and  parallel  to  a  tolerance  of  jhO.OOl  inch. 

•  A  section  of  each  specimen  was  preserved  for  microscopic  observation.  TTiese  virgin 
specimens  were  used  as  a  reference  for  comparison  with  post-test  specimens. 

•  The  specimens  were  carefully  weighed  and  dimensions  recorded.  Each  specimen  was 
instrumented  with  strain-gauged  cantilevers  for  strain  measurements. 

•  Each  specimen  was  installed  into  a  servo-controlled  stiff  testing  machine.  Strain 
rates  for  the  tests  were  low  (10'^  in/in/sec  to  5  x  10*^  in/in/sec)  to  assure  pseudo-static 
conditions.  With  strain  control,  the  complete  stress-strain  curve  (including  the  post-failure 
region)  was  readily  reordered. 

•  At  the  stress  or  strain  of  interest  (linear  elastic  deformation  region,  microcrack 
generation  region,  or,  post-failure  region  on  a  stress-strain  curve),  the  test  was  stopped 
and  specimens  were  unloaded.  Unloading  was  done  slowly  to  avoid  introducing  additional 
microcracks. 

•  To  preserve  the  microstructure  of  the  specimens,  petrographic  epoxy  or  low  melting 
point  rnetal  was  injected.  The  low  melting  point  metal  injection  technique  has  been  used 
very  successfully  for  preserving  microstructures  (Zheng  et  al.  1989  and  Zheng,  1989). 
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With  a  surface  tension  of  400  mN/m,  low  melting  point  metal  such  as  Wood’s  metal 
(melting  point  at  89°C)  will  penetrate  cracks  with  apertures  of  0.08  micrometer,  as  a  liquid 
under  a  pore  pressure  of  10  MPa. 

•  These  specimens  were  sectioned  parallel  and  perpendicular  to  their  axes.  Sections 
for  microscopic  observation  using  scanning  electron  microscopy  and  petrographic 
microscopy  (thin  sections)  were  prepared. 

•  From  the  microscopic  observations,  the  characteristics  of  stress-induced  microstruct¬ 
ures  in  the  specimens  will  recorded.  This  will  include  the  density,  length  distribution, 
aperture  distribution  and  orientation  of  the  microcracks  and  fractures. 

•  These  characteristics  of  the  stress-induced  microstructures  were  statistically  analyzed 
in  association  with  the  stress  conditions  to  which  the  specimen  was  subjected. 

•  Results  of  the  statistical  analysis  of  such  stress-microstructure  relations  provided  the 
basis  for  the  effective  elastic  property  analysis  in  the  theoretical  portion 

of  this  proposal  (Task  2). 

Output:  Statistical  representation  of  the  distribution  of  microcracks  generated  under 
different  stress  conditions.  These  include  the  density,  orientation,  aperture  and  length 
distribution  of  microcracks,  as  functions  of  stress  conditions. 


23  Task  2:  Theoretical  Study  of  Effective  Elastic  Properties  as  Functions 
of  Microfractures 

Objective:  In  this  task,  the  anisotropic  effective  elastic  properties  due  to  stress  induced 
microstructures  was  analytically  derived. 

Technical  Approach:  The  basis  of  this  theoretical  study  is  the  principle  of  energy 
conservation.  The  elastic  strain  energy  derived  using  effective  mechanical 
properties  should  balance  the  elastic  strain  energy  used  in  deformation  of  the  intact  portion 
of  a  fractured  material  and  the  fracture  energy  in  opening  or  closing  of  the  fractures: 

The  derived  anisotropic  elastic  properties  was  validated  qualitatively  using  the  experimental 
results  (stress-strain  curves  and  statistical  distribution  of  microstructures  induced  by 
different  stress  conditions)  from  Task  1  and  artificial  material  by  the  means  of  non-destru¬ 
ctive  (acoustic)  elastic  property  measurements. 

Output:  Anisotropic  effective  elastic  properties  as  functions  of  microcrack  distributions 
(density,  orientation,  length  and  aperture),  summary  of  theoretical  considerations  in 
interrelating  micro-  and  macro-  characteristics. 
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3  Introduction 


Mechanical  properties  of  a  material  change  when  applied  stresses  become  high  enough  to 
induce  microstructural  alterations.  This  is  specially  true  for  geomaterials,  which  are 
generally  heterogeneous  and  are  often  brittle.  Interaction  between  these  microstructural 
alterations  can  ultimately  cause  macroscopic  failure  of  the  materials.  During  the 
occurrence  of  these  microstructural  alterations  (damage),  the  macroscopic  mechanical 
behavior  of  a  material  may  also  change.  Such  evolution  is  usually  represented  as  changes 
in  the  elastic  modulus  of  the  material.  To  appreciate  this,  consider  Figure  1.  This  is  a 
“typical”  axial  deviatoric  stress-strain  curve  generated  during  triaxial  compression.  The 
curve  shows  four  generic  portions  or  segments.  In  segment  A,  initial  large  deformation  due 
to  the  closure  of  microcracks  and  pores  is  characteristic.  The  slope  of  the  stress-strain 
curve  (elastic  modulus)  increases  with  increases  in  the  axial  stress  as  compaction  of  the 
material  proceeds,  reaching  a  constant  value.  The  slope  in  this  linear  portion  (segment  B) 
is  usually  taken  as  the  elastic  modulus  of  the  material.  With  further  elevation  of  the 
deviatoric  stress  (segment  C),  an  alternate  trend  begins,  usually  at  ==  70%  to  80%  of  the 
peak  stress  level.  In  this  regime,  the  stress-strain  curve  manifests  a  decreasing  slope  due 
to  the  generation  and  propagation  of  microcracks,  pore  collapse  and  other  irreversible 
microstructural  changes.  In  segment  D,  generally  regarded  as  the  post-failure  or  strain 
softening  regime,  deformation  increases  with  decreasing  axial  stress.  In  post-peak  regimes 
characteristic  phenomena  are  due  to  the  interaction  of  microcracks  and  large  scale 
microstructural  alterations.  An  understanding  of  the  relationship  between  the  microstruc¬ 
tural  alterations  and  the  macroscopic  behavior  of  a  material  is  important  in  understanding 
mechanisms  of  geomaterial  failure. 

The  behavior  of  rocks  under  compressive 
stresses  has  been  studied  by  vast  numbers  of 
researchers.  For  example,  Wawersik  and 
Fairhurst  (1971)  studied  the  macroscopic 
behavior  of  different  types  of  rocks  in  the 
post-failure  region.  Hallbauer  et  al.  (1973) 
studied  microdamage  density  in  rocks  along  a 
stress-strain  curve  and  showed  that  the  density 
of  microcracks  increases  as  the  applied  stress 
increases.  Fonseka  et  al.  (1985)  studied  micro¬ 
crack  development  under  compressive  stress 
using  scanning  electron  microscopy  and  acous¬ 
tic  emission.  Zheng  et  al.  (1989)  studied 
microcrack  distribution  as  a  function  of  confin¬ 
ing  stress  conditions  using  low  melting  point 
alloy  injection  techniques.  Most  of  these 
studies  were  performed  on  rock.  The  mea¬ 
surements  and  observations  required  destruc¬ 
tion  of  the  material  (except  for  the  technique  used  by  Fonseka  et  al.  1985),  and  the 
information  obtained  is  only  applicable  to  the  stress  state  under  which  the  material 


Figure  1.  Complete  stress-strain  curve 
showing  different  stages  of  rock  behavior. 
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microstructure  was  preserved.  The  research  presented  here  was  initiated  to  observe  and 
record  damage  of  Indiana  limestone  and  concrete  (with  small  aggregate)  at  different  stages 
of  triaxial  compression.  Methods  used  included  concurrent  measurement  of  static  (stress 
and  strain  measurements)  and  dynamic  (ultrasonic  velocity  measurements)  properties  and 
injection  of  a  low  melting  point  metal  alloy. 


4  Experimental  Technique 

Cylindrical  Indiana  limestone  (composition*:  calcium  carbonate  =  97.39%,  magnisium 
carbonate  =  1.2%,  silica  =  0.69%,  alumina  0.44%,  iron  oxide  0.18%  and  porosity  =  10%) 
and  small  aggregate  concrete  specimens,  one  inch  in  diameter  by  nominally  two  inches  long, 
were  endground  parallel  to  within  a  tolerance  of  +.0.001  inch.  These  specimens  were  then 
oven  dried  at  110  °C  for  at  least  48  hours.  Before  testing  ,  each  specimen  was  placed 
between  two  stainless  steel  endcaps  and  jacketed  with  a  Teflon  sleeve,  to  prevent  intrusion 
of  confining  fluid  into  the  specimens.  Each  of  the  endcaps  was  contained  embedded 
ultrasonic  transducers  as  well  as  a  central  axial  hole  and  grooves  on  the  faces  contacting 
the  specimen  for  pore  fluid  access.  Each  specimen  was  then  instrumented  with  strain 
measurement  devices  (strain-gaged  cantilevers,  as  shown  in  Figure  2)  and  placed  in  a 
pressure  vessel. 

To  preserve  microstructure  and  microstructural  alterations  evolving  during  a  test,  a  low 
melting  point  metal  alloy,  Cerrosafe*  was  used  as  the  pore  fluid.  Cerrosafe*  is  composed 
of  42.5%  Bi,  37.7%  Pb,  11.3%  Sn  and  8.5%  Cd,  with  a  melting  point  of  87  °C\  It  has 
been  previously  used  for  porosimetry  studies  (Yadev  et  al.,  1984),  fracture  contact  area 
studies  (Pyrak  et  al.,  1987)  and  microstructure  preservation  (Zheng  et  al,  1987,  Ewy  et  al., 
1989,  Zheng  et  al,  1989,  Suarez  et  al,  1990).  By  injecting  liquid  alloy  into  a  specimen  and 
allowing  it  to  solidify  in  place,  the  microstructure  existing  at  the  instant  of  alloy 
solidification  is  preserved.  TTie  technique  not  only  preserves  microstructure  under  applied 
stress  or  strain  boundary  conditions,  but  also  provides  a  means  of  distinguishing  these 
discontinuities  from  those  generated  during  unloading  and  subsequent  material  handling. 
With  a  surface  tension  of  400  mN/cm^  the  alloy  can  penetrate  fractures  with  apertures  of 
as  little  as  0.08  mm,  with  an  injection  (pore)  pressure  of  10  MPa  (1450  psi). 

Before  testing,  the  temperaiure  of  the  vessel  was  elevated  to  92°C.  A  hydrostatic  stress  of 
2000  psi  was  applied  to  the  specimen.  Next,  pore  pressure  was  elevated  to  1500  psi.  An 
equilibrium  period  was  imposed  to  insure  that  the  alloy  penetrated  every  accessible  void 
in  the  specimen.  The  confining  pressure  was  then  increased  to  the  required  level  for 
testing.  Next,  a  deviatoric  (axial)  stress  was  increased.  This  axial  deviatoric  stress  was 
increased  under  strain  (axial)  control,  at  a  rate  of  10'^  in/in/sec. 


*  Handbook,  Indiana  Limestone  Institute  of  America,  Inc.,  17th  edition,  1986. 
^  Cerro  Alloys  Application  Data,  Cerro  Metal  Products  Co. 
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In  order  to  examine  elastic  property  change  as  a  function  of  microstnictural  evolution 
“continuously”  on  the  same  sample,  concurrent  ultrasonic  velocity  measurements  were 
performed.  Microscopic  observations  only  allow  assessment  of  the  state  of  material  at  the 
moment  of  preservation.  The  endcaps  used  during  measurements  contained  embedded 
longitudinal  (P)  wave  and  shear  (S)  wave  transducers  for  emission  and  receiving  of 
ultrasonic  impulses.  As  testing  proceeded,  the  velocities  of  these  P  and  S  waves  were 
periodically  measured.  Dynamic  Young’s  modulus  and  Poisson’s  ratio  were  calculated  from 
the  velocities  and  stress-strain  curves  were  artificially  reconstructed  from  these  calculated 
elastic  (dynamic)  properties.  The  stress-strain  curves  were  later  compared  with  those 
generated  from  direct  measurement  of  stress  and  strain. 

The  alloy  impregnated  specimens  were  sectioned  parallel  to  their  axes  and  surface  polished 
for  microscopic  observation  and  statistical  studies.  Microcrack  density,  average  orientation, 
length  and  aperture  were  recorded  in  the  same  manner  as  described  by  Zheng  et  al;  (1989). 
The  statistical  results  were  then  used  in  effective  property  measurements  and  evaluation. 


5  RESULTS 


5.1  “Static”  Stress-Strain  Curves 

Figures  3  to  7  are  typical  stress-strain  curves  for  samples  tested  under  different  confining 
stress  conditions.  Figures  3,  4  and  5  show  the  stress-strain  behavior  of  Indiana  limestone 
under  confining  stresses  of  0,  2500  and  5000  psi,  respectively.  Figures  6  and  7  show  the 
stress-strain  behavior  of  concrete  specimens  subjected  to  confining  stresses  of  0  and  5000 
psi. 


5.2  “Dynamic”  Stress-Strain  Curves 

Figures  8  to  12  are  analogous  stress-strain  curves  artificially  constructed  fi-om  moduli 
determined  fi’om  the  concurrent  ultrasonic  measurements.  These  figures  correspond 
directly  with  Figures  3  through  7,  performed  at  a  quasi-static  strain  rate.  Detailed  data  is 
provided  in  Appendix  A.  The  “dynamic”  stress-strain  curves  were  constructed  by  using 
dynamic  modulus  as  the  secant  modulus  of  a  stress-strain  curve  at  corresponding  stress 
levels^.  The  dynamic  modulus  and  Poisson’s  ratio  were  calculated  using  the  following 
equations  (Jaeger  and  Cook,  1976): 


^  Damage  mechanics  considerations  frequently  adopt  secant  modulus  because  of  the 
inconsistency  of  tangential  moduli  in  post  peak  regions  (negative). 
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Figure  3.  Quasi-static  stress-strain  plot  for  an  Indiana  limestone  sample,  oriented  vertically 
and  subjected  to  zero  confining  stress.  The  periodic  levels  of  constant  stress  difference 
indicate  where  dynamic  measurements  were  performed. 
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Figure  4.  Quasi-static  stress-strain  plot  for  an  Indiana  limestone  sample,  oriented  vertically 
and  subjected  to  a  confining  stress  of  2500  psi.  The  periodic  levels  of  constant  stress  differ¬ 
ence  indicate  where  dynamic  measurements  were  performed. 
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Figure  5.  Quasi-static  stress-strain  plot  for  an  Indiana  limestone  sample,  oriented  vertically 
and  subjected  to  a  confining  stress  of  5000  psi.  The  periodic  levels  of  constant  stress  differ¬ 
ence  indicate  where  dynamic  measurements  were  performed. 
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Figure  6.  Quasi-static  stress-strain  plot  for  a  concrete  sample,  subjected  to  zero  confining 
stress.  The  periodic  levels  of  constant  stress  difference  indicate  where  dynamic  measure¬ 
ments  were  performed. 


cn 

Q. 


UJ 

CJ 


UJ 

cr 

ui 

u. 


a 

cn 

cn 

UI 

QC 

I— 

cn 


I  TerraTek  |  (RADIAL)  STRAINS.  X  (AXIAL) 

Figure  7.  Quasi-static  stress-strain  plot  for  a  concrete  sample,  subjected  to  a  confining  stress 
of  5000  psi.  The  periodic  levels  of  constant  stress  difference  indicate  where  dynamic 
measurements  were  performed. 
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Axial  Strain  (%) 


Figure  8.  “Dynamic”  stress-strain  curve  plot  for  an  Indiana  limestone  sample,  oriented 
vertically  and  subjected  to  zero  confining  stress.  The  plot  was  constructed  from  dynamic 
modulus  calculated  from  the  ultrasonic  wave  velocities. 


Axial  Strain  (%) 

Figure  9.  “Dynamic”  stress-strain  curve  plot  for  an  Indiana  limestone  sample,  oriented 
vertically  and  subjected  to  a  confining  stress  of  2500  psi.  The  plot  was  constructed  from 
dynamic  modulus  calculated  from  the  ultrasonic  wave  velocities. 
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Axial  Strain  (%) 


Figure  10.  “Dynamic”  stress-strain  curve  plot  for  an  Indiana  limestone  sample,  oriented 
vertically  and  subjected  to  a  confining  stress  of  5000  psi.  The  plot  was  constructed  from 
dynamic  modulus  calculated  from  the  ultrasonic  wave  velocities. 
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Figure  11.  “Dynamic”  stress-strain  curve  plot  for  a  concrete  sample,  oriented  vertically  and 
subjected  to  a  zero  confining  stress.  The  plot  was  constructed  from  dynamic  modulus 
calculated  from  the  ultrasonic  wave  velocities. 
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Axial  Strain  (%) 

Figure  12.  “Dynamic”  stress-strain  curve  plot  for  a  concrete  sample,  oriented  vertically  and 
subjected  to  a  confining  stress  of  5000  psi.  TTie  plot  was  constructed  from  dynamic 
modulus  calculated  from  the  ultrasonic  wave  velocities. 


S3  Comparison  of  Static  and  Reconstructed  Dynamic  Stress-Strain 
Curves 

In  general,  the  reconstructed  stress-strain  curves  show  smaller  axial  strains.  This  is 
consistent  with  the  tendency  for  dynamic  moduli  to  be  higher  than  static  moduli.  This  is 
not  the  issue  of  import  here.  This  phenomena  is  well  established  in  the  literature  (King, 
1980,  for  example).  The  issues  of  major  importance  is  a  qualitative  determination  of 
whether  dynamic  measurements  are: 

(i)  An  accurate  reflection  of  the  relative  number  of  microcracks  present,  and/or, 

(ii)  Indicate  catastrophic/permanent  deformation. 

As  acoustical  testing  is  commonly  used  to  assess  the  presence  of  flaws,  this  is  a  particularly 
serious  exercise.  Initial  observations  tend  to  suggest  that  the  shapes  of  the  static  and 
dynamic  curves  are  similar.  The  indicators  are  qualitative  only  in  this  Phase  I  effort. 


5.4  Photomicrograph 

Figures  13a  and  13b  through  Figures  17a  and  17b  are  micrographic  representations  of  the 
tests  illustrated  previously.  In  each  case  (a)  is  a  photograph  at  a  magnification  of 
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approximately  12  times  (the  width  of  the  photograph  shown  is  ore  inch)  and  (b)  is  a  trace 
of  dominant  microcrack  features.  A  further  observation  relates  to  the  differing  microcrack 
spectrums,  particularly  different  averages  orientation  between  the  limestone  and  the 
concrete.  Table  1  details  simple  statistical  processing  of  microcrack  data  from  these 
samples,  as  a  function  of  hydrostatic  (and  deviatoric)  stress  conditions. 


5.5  Microcrack  Observations 

From  Figures  13  through  17,  it  can  be  seen  that  most  of  the  microcracks  are  extensile  and 
are  parallel  or  subparallel  (small  angle)  to  the  direction  of  the  maximum  applied 
compressive  stress.  Table  1  shows  that  the  average  absolute  angle  between  microcrack 
planes  and  the  direction  of  the  maximum  applied  compressive  stress  is  less  than  15°.  The 
standard  deviation  for  the  angle  with  respect  to  the  sample  axis  for  all  of  the  tested  samples 
is  ~  10°.  In  the  limestone  samples,  the  photomicrographs  show  that  most  microcracks 
propagate  linearly  except  near  their  tips  where  they  turn  into  slip  planes  of  calcite  crystals, 
just  before  they  terminate. 

Microcracks  in  the  limestone  appear  to  have  been  generated  by  more  than  one  mechanism. 
However,  these  different  mechanisms  all  relate  to  local  tensile  stress.  Many  microcracks 
were  generated  due  to  bending  of  long,  beam-like  grains  (with  their  long  axes  parallel  to 
bedding).  In  addition,  contact  point  loading  generated  tensile  stress  concentrations  which 
also  generated  extensile  features.  Tensile  stress  concentrations  around  an  open  pore  space 
or  a  soft  inclusion  can  also  generate  microcracks.  Microcracks  were  found  to  have  been 
generated  from  grain  boundaries,  especially  those  parallel  to  the  direction  of  the  maximum 
compression.  In  the  concrete  samples,  most  of  the  microcracks  were  generated  from  pore 
spaces  and  grain  (aggregate)  boundaries. 

The  density  of  the  microcracks  appears  to  increase  as  the  confining  stress  is  increased.  The 
length  distribution  of  microcracks  is  strongly  dependent  on  the  confining  stress  applied. 
At  zero  and  low  confining  stress,  microcracks  apparently  propagated  farther  than  those 
generated  under  higher  confining  stress,  before  they  were  arrested.  Interaction  between 
microcracks  occurs  more  readily  at  low  confining  stress  because  the  microcracks  are  longer. 
This  resulted  in  sample  failure  at  relatively  lower  axial  stresses.  The  decrease  in  microcrack 
length  as  the  confining  stress  was  increased  is  due  to  the  stress  intensity  at  crack  tips.  As 
indicated  above,  microcracks  are  generated  by  local  tensile  stresses  which  are  usually 
concentrated  at  the  origin  of  the  microcracks.  Propagation  of  these  microcracks  is  driven 
by  this  local  tensile  stress.  The  confining  stress  (global  compression)  generates  a  negative 
stress  intensity  component  counteracting  the  stress  intensity  factor  due  to  the  local  tensile 
stress.  As  a  result  in  these  tests,  confining  stress  reduces  the  stress  intensity  factor  that 
propagates  the  microcracks.  Hence,  the  average  length  of  microcracks  was  reduced  with 
an  increase  in  confining  stress. 

In  general,  crack  aperture  (width)  appears  to  be  proportional  to  crack  length,  since  the 
aperture  of  a  crack  is  also  proportional  to  stress  intensity  at  crack  tips.  As  can  be  seen 
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Figure  13b.  Trace  of  dominant  microcracks  shown  on  the  photomicrograph  in  Figure  13a. 
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I  Figure  14a.  Photomicrograph  of  an  axial  cross  section  Indiana  limestone  sample,  subject¬ 

ed^  to  a  confining  stress  of  2500  psi  (post-peak).  The  width  of  the  photomicrograph  is 
the  diameter  of  the  sample  and  is  one  inch  in  dimension. 
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I  I  mure  15a.  Photomicrograph  of  an  axial  cro.s.s  section  Indiana  limestone  sample,  subject- 

I  eef  to  a  confining  stress' of  .5000  psi  (post-peak).  The  width  of  the  photomicrograph  is 

I  hi  diameter  of  the  sample  and  is  one  inch  in  dimension. 
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1  igure  16a.  Photomicrograph  of  an  axial  cross  section  of  concrete  sample,  subjected  to 
a  zero  confining  stress  (post-peak).  The  width  of  the  pliotomicrograph  is  the  diameter 
Ilf  the  sample  and  is  one  inch  in  dimension. 
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F  igure  17a.  Pliotoniicrograph  of  an  axial  cross  section  of  concrete  sample,  subjected  to 
a  confining  stress  of  .‘'OOO  psi  (post-peak).  The  width  of  the  photomicntgraph  is  the 
diameter  of  the  sample  and  is  one  inch  in  dimension. 
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Figure  17b.  Trace  of  dominant  microcracks  shown  on  the  photomicrograph  in  Figure  17a. 

I 
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from  Table  1,  the  crack  aperture  decreased  with  increasing  confining  stress.  It  appears  that 
the  average  aspect  ratio  of  microcracks  (width  to  length  ratio)  in  the  limestone  varies 
substantially  but  it  was  very  consistent  in  concrete.  This  is  a  characteristic  difference 
between  the  limestone  and  the  concrete  samples. 

Microcrack  interaction  in  the  limestone  was  observed  to  occur  in  localized  pattern.  Under 
unconfined  compression  (Figures  3,  8  and  13),  microcracks  tended  to  connect  with  each 
other,  form  a  longer  extensile  fracture  which  usually  caused  axial  splitting.  Under  triaxial 
compression,  where  shear  failure  is  in  general  the  typical  failure  mode,  it  was  observed  that 
shear  fractures  were  formed  by  many  small  extensile  microcracks  parallel  to  and  offset  from 
each  other  (Figures  5,  10  and  15).  In  concrete,  however,  no  macro  shear  surfaces  were 
observed  and  all  samples  failed  in  axial  splitting.  Most  of  the  fractures  were  composed  of 
a  few  smaller  cracks  but  not  as  many  as  in  the  limestone. 


5.6  “Dynamic”  Stress-Strain  Curves 

The  “dynamic”  stress-strain  curves,  obtained  from  concurrent  ultrasonic  velocity 
measurements,  were  rather  similar  to  those  generated  from  conventional  static  stress-strain 
monitoring.  It  appears  that  these  dynamic  curves  do  not  show  the  initial  compaction 
evident  in  the  static  stress-strain  curves.  However,  the  “dynamic”  stress-strain  curves  do 
show  the  generation  and  increase  in  density  of  microcracks.  The  measured  changes  in 
acoustic  velocities  are  significant  because  they  reflect  the  microstructural  changes  in  the 
rock  or  the  concrete.  In  combination  with  acoustic  emission  methods  (recording  acoustic 
events  caused  by  generation  of  microcracks  in  a  testing  sample),  it  appears  feasible  that 
changes  in  the  dynamic  mechanical  properties  may  be  used  for  determination  of  microcrack 
density. 
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SBIR  Microcrack  Population 


6  THEORETICAL  CONSIDERATIONS  ON  MECHANI¬ 
CAL  PROPERTY-MICROCRACK  RELATIONSHIP 


6.1  Background 

The  previous  section  on  experimental  observations  provides  qualitative  representation  of 
damage  evolution.  Beyond  empirical  observations,  it  is  desirable  to  quanti^  the  influence 
of  “microcracks”  in  bulk  material  behavior.  This  is  done  in  the  following  section, 
presuming  elastic  relationships  and  various  presumptions  on  axes  of  symmetry. 


6.2  CONSTITUTIVE  EQUATIONS 


An  orthotropic,  homogeneous  and  linear  elastic  material  has  different  mechanical 
properties  in  three  orthogonal  directions.  The  constitutive  relationship  for  such  a  material 
can  be  expressed  in  a  1,  2,  3  Cartesian  coordinate  system  as  follows: 
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where  and  Ojj  are  the  strain  and  stress  components  on  the  i  plane,  in  the  j  direction;  Vij 
is  the  deformation  in  i  the  direction  induced  by  stress  in  the  j  direction  due  to  Poisson’s 
effect;  Ej  and  G;  are  Young’s  modulus  and  shear  modulus  in  the  i  direction. 


Equation  (1)  can  be  expressed  in  matrix  form: 
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The  inverse  of  the  compliance  matrix  (A),  often  regarded  as  the  stiffness  matrix  (C)  can 
be  expressed  as: 


(o)=(A)->(€)=(C)(e)= 
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This  matrix  exists  if  and  only  if  the  compliance  matrix  is  symmetric.  As  a  result,  the 
following  relationships  exist: 
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Matrices  (C„)  and  (C,)  are  the  inverse  matrices  of  (A„)  and  (A,),  respectively: 
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6.3  EQUILIBRIUM  EQUATIONS'* 

The  equilibr.uin  equations  for  an  elastic  material  in  the  absence  of  body  forces  can  be 
written  as  foilows: 


*For  more  information,  refer  to  Jaeger  and  Cook,  1976. 


29 


u^=B^XjCosa^  +Ct) 

In  equation  (8),  Bj  is  a  constant,  Qfj  is  the  angle  of  incidence  of  an  elastic  wave  with  respect 
to  the  normal  of  plane  i  (direction  i),  and  C  is  the  velocity  of  the  elastic  wave. 
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(9) 


6.4  ELASTIC  WAVE  MOTION  EQUATION 

By  bringing  (3)  and  (7)  into  (6),  one  can  obtain  the  motion  equation  for  an  elastic  wave: 

^21  ^21  '^^72  ^72*^73  ^23~P^^^2 
^31  ^31^^32  ^32^^33  ^33“P^^^3 
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Dji  =CjjCOS^a  j  +G3COs^a2+G2COS^a3 
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1)33  =C33Cos^a2 +GjCOS^a2 +G2COS^O  j 
Dj2=(Cj2+G3)cosajCosa2 
D23 =(C23 +Gj)cosa2Cosa3 
^31  ^(^31  '*^G2)^osa^cosa  j 

The  wave  motion  equations  (9)  can  be  written  in  matrix  form: 
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where  A  =  pC^ 


A  set  of  non-zero  solutions  for  (10)  or  (11)  exists  if  and  only  if  the  determinant  of  the 
coefficient  matrix  of  equation  (11)  equals  zero,  i.e., 
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The  values  of  k  which  satisfy  equations  (12)  or  (13)  are  eigenvalues  of  equation  (11).  The 
velocities  C  of  elastic  waves  can  be  obtained  from  C=V  Ayp ,  once  the  eigenvalues  k  are 
obtained. 


6.5  ELASTIC  RELATIONS  FOR  AN  AXISYMMETRICALLY  MICROCR- 

ACKED  MATERIAL 

Relations  outlined  in  previous  sections  are  well  known.  The  following  derivations 
incorporate  microcracking.  Restricting  consideration  to  microcracks  generated  in  a 
homogeneous,  isotropic  and  elastic  material  under  triaxial  compression  (aj  >  a^),  the 

elastic  properties  of  this  cracked  (damaged)  material  can  be  considered  to  be  axisymmetric 
and  the  elastic  modulus  in  the  axial  direction  is  presumed  not  to  be  changed.  The  elastic 
properties  of  this  material  are: 

£j=£  (known,  same  as  for  the  original  uncracked  (undamaged)  material),  £2“^3 
V2j=V3|=v(known,  sameastheuncracked(undamaged)material),  v,2=Vj3,  V23=V32  (1^) 
Gj,  G2=G3 

Using  (14)  in  the  constitutive  equations  (2a)  and  (2b),  and  (5a)  and  (5b);  one  can  obtain: 
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6.6  ELASTIC  CONSTRAINT  (ENERGY  CONSERVATION)  FOR  AN 
ISOTROPIC  MATERIAL  WITH  AN  ANISOTROPIC  EXTENSILE 
MICROCRACK  DISTRIBUTION 

The  principles  of  superposition  and  energy  conservation  can  be  used  to  derive  the  elastic 
constants  for  an  isotropic  material  containing  preferentially  distributed  (oriented)  extensile 
cracks  (Figure  18),  as  defined  in  section  4.4. 

Using  energy  conservation  principles: 

VJ  =  U^*Nu^  (16) 

where  U'  is  the  total  strain  energy  in  the  cracked  (damaged)  material,  U,  is  the  total 
energy  for  the  uncracked  material,  N  is  the  total  number  of  cracks  and  u^  is  the  strain 
energy  of  a  single  crack.  U'  and  are  functions  of  the  applied  stress  and  elastic 
properties  of  the  material: 
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Figure  18.  Schematic  representation  of  the  principles  of  superposition  applied  to  a  material 
containing  microcracks. 
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and 

U=^oe{E,\)  (17b) 

where  V  is  the  volume  of  the  material  and  all  the  variables  in  parentheses  are  elastic 
properties  of  a  cracked  material  (17a)  and  the  original  uncracked  material  (17b). 

Presume  that  typical  microcrack  morphology  is  penny  shaped.  The  strain  energy  around 
a  penny-shaped  crack  (Figure  19)  can  be  expressed  as  the  energy  per  unit  area  of  the  crack: 

j ,  (18) 

E  ^  E  ^  E 

where  J  is  the  energy  released  per  unit  crack  surface  area  generated,  K„  K,„  and  K,„  are 
the  stress  intensity  factors  for  Mode  I,  II  and  III  crack  surface  displacements. 
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Using  (17)  and  (18)  into  (16)  one  obtains: 


34 


(19) 
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where  A  is  the  surface  area  of  a  single  crack. 


6.7  STRAIN  ENERGY  FOR  A  PENNY-SHAPED  CRACK 

The  stress  intensity  factors  for  a  penny-shaped  crack  under  uniform  tensile  and  shear  stress 
(Figure  19)  are  as  follows: 
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where  a  is  the  radius  of  the  penny-shaped  ciack,  r^^/lpr  Tjj)  is  the  uniform  shear  stress 
acting  on  planes  2  (or  3)  in  the  direction  of  1,  ctq  is  the  uniform  tension  orthorganal  to  the 
crack  plane,  T23  is  the  uniform  shear  stress  on  the  crack  plane  and  9  is  the  on  the  crack 
plane  angle  from  direction  2  to  the  point  of  consideration. 
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Figure  19.  Penny-shaped  crack  in  an  infinite  body  under  uniform  tension  and  shear. 
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The  strain  energy  around  this  single  crack  is: 
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These  energy  relationships  will  be  used  in  subsequent  sections  for  delineating  effective 
modulus. 


Using  (20)  in  (21),  one  obtains: 
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6.8  ELASTIC  (P  AND  S)  WAVE  PROPAGATION 

During  compressional  (P)  wave  prc  pagotion  through  an  elastic  medium,  the  medium  is 
effectively  subjected  to  uniaxial  elasii"  scrain,  that  is,  if  the  P-wave  is  parallel  to  direction 
i,  then, 

€..i«s0,  €„  =  e^^  =  0,;Vi,  kT^i 

In  the  following  subsections,  P  and  S  wave  propagation  in  a  material  with  axisymmetrically 
distributed  cracks  (generated  from  triaxial  compression)  are  evaluated.  Wave  propagation 
in  different  principal  directions  is  considered. 
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The  Dormal  and  shear  stresses  applied  on  a  crack  whose  normal  is  oriented  at  (j)  from 
direction  2  are: 


k*=o 


(29) 


Assuming  that  N  cracks  are  homogeneously  distributed  with  angle  (f>  varying  from  0  to  360° 
(axisymmetric),  the  total  strain  energy  for  the  microcracks  can  be  expressed  as: 


The  energy  conservation  equation  (19)  yields: 
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Using  equations  (27),  (28),  (29),  and  (22)  in  (31),  gives: 
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These  equations  will  be  used  later  to  solve  for  effective  moduli  and  Poisson’s  ratios. 


6.8.2  Shear  (S)  Wave  Traveling  in  Direction  1: 

(Particle  Motion  in  a  Plane  with  its  Normal  Parallel  to  Direction  1) 

Assuming  that  a  shear  (S)  wave  (with  magnitude  t)  travels  in  direction  1  and  vibrates  with 
an  angle  a  to  direction  2,  the  stresses  on  a  crack  at  an  angle  to  direction  2  are: 

The  energy  conservation  principle  (19)  gives: 

^  32Na^(l-«-v)(2-t-v^)^2 
2Gj~2G^  3KB(2-v)  ^ 

1  _  1  ,32x(l^v)(2^v^) 

2C?j  2G  3£(2-v) 

where 

No  ^ 

X= - =  crack  density  (dimensionless) 


6.8.3  P  and  S  Wave  Traveling  Orthogonal  to  Direction  1 

The  following  example  presumes  a  wave  propagating  in  direction  2.  Derivations  for  other 
“radial”  directions  are  analogous  because  of  symmetry  conditions.  When  a  P-wave  (with 
magnitude  a^)  is  propagating  in  direction  2,  the  material  is  subjected  to  uniaxial  strain  in 
direction  2. 
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When  a  shear  (S)  wave  is  traveling  (with  magnitude  t)  in  direction  2  and  vibrates  with  an 
angle  a  to  direction  1,  the  stress  components  Tij,  213,  and  T23  are 


Tj2=xcosa 

Tj3=xsina 

t,3=0 


Stresses  acting  on  a  fracture  plane,  refer  to  Figure  20,  are: 
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The  energy  conservation  equation  becomes: 
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Using  equations  (38),  (39)  and  (40)  in  (41)  and  (42),  one  obtains: 
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Figure  20.  Stresses  acting  on  a  penny-shaped  crack  whose  normal  is  orthogonal  to  direction 
1  and  at  an  angle  (f>  to  direction  2. 


40 


Because  the  incidence  angle  a  varies  from  0  to  tt,  (44)  becomes  (after  integration  with 
respect  to  a): 
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Solving  (33),  (36),  (43)  and  (45)  simultaneously  gives  the  effective  mechanical  properties 
of  a  multiple-cracked  material: 
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Equations  (46)  to  (50)  are  expressions  for  effective  moduli  and  Poisson’s  ratio  for  an 
axisymmetrically  microcracked  material,  as  functions  of  the  original  (undamaged)  material 
properties  and  crack  density.  In  other  words,  if  the  elastic  properties  of  an  initially 
homogeneous  material  and  the  density  of  microcracks,  which  are  generated  by  axial 
compression  are  known,  the  effective  moduli  and  Poisson’s  ratio  in  principal  directions  can 
be  obtained  from  these  equations. 

Relationships  given  in  equations  (46)  through  (50)  also  provide  a  means  of  detecting  micro¬ 
crack  (damage)  density  when  the  effective  properties  can  be  measured  either  statically  or 
dynamically.  In  Phase  II  of  this  research,  relationships  between  dynamic  properties  and 
microcrack  density  will  be  established  and  experimentally  validated.  Non-destructive 
(ultrasonic,  transmitted  waves)  methodology  will  be  proposed  as  a  means  for  detecting  and 
quantifying  damage  (microcrack  density)  i  vel. 


7  Conclusions  and  Observations 

Detailed  microstructiiral  alterations  (damage)  in  Indiana  limestone  and  one  type  of 
concrete  under  different  loading  conditions  have  been  observed  using  liquid  metal  alloy 
preservation  techniques.  The  basic  observations/conclusions  are  as  follows; 

•  As  has  been  cited  elsewhere,  the  basic  mechanism  for  extensile  microcracks  to 
initiate  in  a  material  under  compression  is  the  localized  tensile  stresses  induced  by 
heterogeneities.  In  Indiana  limestone,  the  mechanisms  for  the  generation  of  these  localized 
tensile  stresses  include  bending  of  long  grains,  point  contact  between  adjacent  grains, 
mismatching  elastic  properties,  boundary  discontinuities,  pore  space  and  soft  inclusion 
induced  tension.  In  the  concrete  samples,  however,  most  of  the  microcracks  are  generated 
at  grain  (aggregate)  boundaries. 

•  Once  the  microcracks  are  generated,  sometimes  at  a  large  inclination  angle  with 
respect  to  the  maximum  compression,  they  generally  propagate  in  the  direction  of  the 
maximum  compression.  The  propagation  of  these  microcracks  is  controlled  initially  by  the 
amount  of  local  tensile  stress  responsible  for  the  generation  of  the  microcracks  and  the 
overall  confining  stress,  before  any  interaction  takes  place.  The  generation  and  propagation 
of  microcracks  are  speculated  to  also  change  the  local  stress  distribution  which  in  turn 
affects  the  generation  of  other  microcracks. 

•  Preliminary  statistical  results  for  both  the  concrete  and  limestone  show  that  the 
average  length  of  the  microcracks  decreases  with  an  increase  in  confining  stress.  Average 
aperture  of  the  microcracks  also  decreases  with  an  increase  in  confining  stress.  The  spread 
of  both  length  and  aperture  distribution  decreases  with  an  increase  in  confining  stress. 
However,  the  aspect  ratio,  the  ratio  between  the  aperture  of  a  crack  and  its  length,  remains 
relatively  constant  for  concrete  and  fluctuates  for  Indiana  limestone.  Most  of  the 
microcracks  are  parallel  or  subparallel  to  the  direction  of  the  maximum  compressive  stress. 
The  average  orientations  of  microcracks  in  all  of  the  samples  tested  are  within  15  degrees 
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deviation  from  the  maximum  compressive  direction.  The  density  of  microcracks,  the 
number  of  microcracks  per  unit  observed  area,  increases  with  an  increase  in  confining 
stress.  This  may  be  due  to  the  higher  deviatoric  stress  which  was  required  to  observe  any 
macroscopic  alteration  of  material  properties.  The  relationship  between  microcrack 
distribution  and  stresses  applied  could  not  yet  be  quantified  in  this  phase. 

•  The  macroscopic  failure  of  Indiana  limestone  samples  under  unconfined  conditions 
is  usually  in  a^ial  splitting,  and  under  high  confining  in  form  of  shear  bands.  The  axial 
splitting  of  specimens  are  usually  the  result  of  one  or  several  extensile  cracks  propagating 
all  the  way  to  *he  ends  of  the  samples.  The  shear  bands  are  observed  to  be  formed  from 
many  short  extensile  microcracks  which  are  subparallel  to  the  direction  of  the  maximum 
compression  and  with  their  centers  slightly  offset  from  each  other.  The  overall  effect  of 
these  extensile  microcracks  is  a  shear  band(s)  (see  Figure  14b). 

•  Based  on  the  results  of  a  limited  number  of  tests,  it  was  not  possible  to  conclusively 
determine  the  relationship  between  dilatation  of  samples,  microcracks  and  pore  collapse. 
However,  it  is  speculated  that  in  the  post-peak  region  or  stress-strain  curves,  the  volumetric 
strain  of  the  tested  samples  is  a  combination  of  microcrack  dilatation  (volume  increase)  and 
pore  collapse  (volume  reduction). 

•  An  increase  in  microcrack  density  decreases  static  and  dynamic  elastic  (secant) 
moduli  (qualitatively).  This  is  deduced  from  continuous  measurement  of  ultrasonic 
velocities.  The  relationship  between  microcrack  distribution  and  stresses  applied  could  be 
quantified  in  future  effort.  “Dynamic”  stress-strain  curves  reconstructed  from  ultrasonic 
velocities  showed  similar  behavior  to  the  concurrently  generated  static  stress-strain  curves. 
This  suggests  that  if  a  relationship  between  the  “dynainic”  and  static  stress-strain  curves  can 
be  established,  non-destructive  method  such  as  ultrasonic  detection  can  be  used  to 
quantitatively  determine  and  detect  the  static  behavior  of  a  geomaterial  and  the  damage 
level. 


•  Theoretical  derivations  of  effective  elastic  properties  of  an  isotropic  material  contai¬ 
ning  axisymmetrically  distributed  microcracks,  using  energy  conservation  principles,  showed 
that  moduli  of  a  material  decrease  with  an  increase  in  damage  density.  This  is  verified  by 
independent  concurrent  triaxial  compression  tests,  ultrasonic  velocity  measurements  and 
microstructural  observations.  These  results  show  that  if  a  one-to-one  relationship  can  be 
established  between  the  theoretically  derived  formulations  and  microcrack  distribution,  it 
could  be  possible  to  extrapolate  microcrack  density  at  any  stage  of  loading.  This  requires 
a  calibration  of  the  theoretical  formulations  given  the  stress  conditions  and  quantification 
of  microcrack  distribution  under  the  same  conditions.  In  a  Phase  II  proposal  for  this 
research  submitted  under  separate  cover,  this  calibration  and  microcrack  density 
extrapolation  is  proposed. 
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In-situ  stress  orientation  and  magnitude  (the  principal  stress  tensor)  is  important 
information  for  designing  excavations  in  rock.  For  example  in  petroleum  industry,  in-situ 
stresses  influence  the  orientation  of  deviated  wellbores  to  increase  production  (by  means 
of  stimulation)  and  to  avoid  stability  problems.  When  a  borehole  is  drilled,  the  tangential 
stress  varies  angularly,  due  to  stress  concentrations  caused  by  removal  of  the  material 
previously  in  the  borehole.  The  variation  of  the  stress  concentrations  is  dependent  on  the 
far-field  (virgin)  stress  orientations  and  magnitudes.  When  the  orientations  of  the  near 
borehole  stress  concentrations  can  be  determined,  it  is  possible  to  determine  the  orientation 
of  the  principal  stresses.  Certain  techniques  such  as  micro-hydraulic  fracturing,  televiewer 
surveys  and  mapping  of  borehole  breakouts  have  been  used  to  deduce  magnitudes  and/or 
orientations.  Each  suffers  some  specific  drawbacks,  either  requiring  prolonged  testing  in 
opeiilioie  or  not  providing  information  which  can  be  definitively  used  to  define  the  stress 
state.  Elastic  wave  propagation  shows  that  the  velocities  of  compressional  and  shear  waves 
travelling  through  a  material  depend  on  the  stress  state  (both  axial  and  confining  stresses). 
In  fact,  velocity  ratios  have  been  used  in  conjunction  with  crude  assumptions  of  no  lateral 
deformation  to  characterize  variation  in  a  pseudo-stress  field.  It  is  proposed  that  if  the 
variation  of  acoustic  properties  circumferentially  around  a  borehole  can  be  quantified  that 
adaptations  of  currently  existing  wireline  log  interpretation  procedures  may  afford  the 
possibility  of  determining  directionality  of  stresses  and  possibly  the  magnitudes. 

Experiments  have  been  performed  on  rock  samples,  to  assess  variation  in  acoustic 
characteristics  as  a  function  of  applied  stress  magnitude.  Cylindrical  Indiana  limestone 
samples  were  oven  dried  and  instrumented  to  measure  deformation  as  a  function  of  applied 
stress.  After  instrumentation,  the  samples  were  installed  in  a  pressure  vessel  and  subjected 
to  varying  levels  of  confining  pressure.  After  equilibration,  the  samples  were  subjected  to 
triaxial  compression  at  a  rate  of  10"'  in/in/s.  Ultrasonic  impulses  were  periodically 
transmitted  through  the  samples.  Compressional  and  shear  wave  velocities  were  measured. 
From  velocity  ratios,  dynamic  mechanical  properties  (Young’s  modulus  and  Poisson’s  ratio) 
were  calculated  and  from  the  variation  of  these  properties  with  deviatoric  stress,  artificial 
“dynamic”  stress  strain  curves  were  reconstructed.  These  artificial  curves  showed  similar 
behavior  to  static  stress-strain  curves  determined  concurrently  in  the  testing.  The  dynamic 
modulus  of  the  sample  increased  with  an  increase  in  the  applied  confining  stress  and  with 
an  increase  in  the  deviatoric  axial  stress.  Near  the  peak  axial  stress,  dynamic  modulus 
started  to  decrease,  indicating  generation  of  microfractures  in  the  samples.  Measuring 
acoustic  velocities  at  different  azimuth,  this  dependence  of  dynamic  modulus  of  rock  on 
applied  stresses  can  be  used  as  a  means  of  determining  in-situ  stress  orientation  and 
possibly  magnitude. 


